Critical heat flux around strongly-heated nanoparticles 
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We study heat transfer from a heated nanoparticle into surrounding fluid, using molecular dy- 
namics simulations. We show that the fluid next to the nanoparticle can be heated well above its 
boiling point without a phase change. Under increasing nanoparticle temperature, the heat flux 
saturates which is in sharp contrast with the case of flat interfaces, where a critical heat flux is 
observed followed by development of a vapor layer and heat flux drop. These differences in heat 
transfer are explained by the curvature induced pressure close to the nanoparticle, which inhibits 
boiling. When the nanoparticle temperature is much larger than the critical fluid temperature, a 
very large temperature gradient develops resulting in close to ambient temperature just radius away 
from the particle surface 

PACS numbers: (68.08.De Liquid-solid interface structure: measurements and simulations; 44.35. +c Heat 
flow in multiphase system;65.80.+n Thermal properties of small particles, nanocrystals, and nanotubes ) 
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Sub-micron scale heat transfer is attracting a growing 
interest [l] , motivated by both fundamental and techno- 
logical issues. The fast emergence of this field is, to a 
large extent, associated with the development of micro 
and nano technologies. In some cases, thermal transfer 
is part of the system function (e.g. the use of nanonuids 
for heat transport or of multilayered materials for ther- 
mal insulation). In other cases, the enhancement of heat 
transfer is a key to a proper operation of the microsys- 
tem (e.g. microprocessors) and involves the integration 
on ever smaller scales of devices such as micro heat pipes. 
Although these systems are of micrometer size, the re- 
gions that limit heat transfer -interfaces, constrictions - 
are often characterized by even smaller lengths, bring- 
ing heat transfer issues into the domain of nanosciences. 
Recent interest in heat transport around nanoparticles 
has arisen in part from the particular properties of the 
so called "nanonuids" @, i.e. colloidal suspensions 
of solid nanoparticles, which exhibit improved thermal 
transport properties. On the fundamental side, a number 
of laser heating studies were performed demonstrating 
even melting of metal nanoparticles without macroscopic 
boiling of the embedding liquid. (4j,|5[. The physics of this 
phenomenon involves a complex interplay between boil- 
ing, heat transfer, and particle-fluid interactions (wet- 
ting), and is still poorly understood. 
In this letter, we use molecular dynamics simulation 
(MD) to study heat transfer around a nanoparticle sur- 
rounded by a volatile fluid. We characterize heat transfer 
in situations where the nanoparticle is heated above the 
fluid boiling point and/or critical temperature, and com- 
pare this situation to the case of flat interfaces. We show 
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that the fluid around nanoparticle can sustain large heat 
fluxes, well above the critical heat flux of the bare fluid 
on a flat interface. Accompanying the large heat flux are 
extreme temperature gradients allowing for localization 
of the hot liquid to volumes comparable with nanoparti- 
cle size. 

MD has been already applied successfully to characterize 
heat transfer across a nanoparticle/fluid interface, in the 
absence of fluid phase change d, 0, [1]. This technique 
has the advantage to give local detailed information on 
heat transfer and on the structure of the fluid close to 
the nanoparticle interface as well. Furthermore, the sys- 
tem at hand is perfectly controlled, which allows to pin- 
point unambiguously the relevant physical mechanisms 
at work. 

The model we simulate consists of a solid nanoparticle 
made of 555 atoms immersed in a fluid of 23000 atoms. 
All atoms interact through a Lennard- Jones potential 
V a p{r) = 4e((er/r) 12 — c a p(<r/r) 6 ) where a, (3 refers to 
solid or liquid atoms. The potential has a cut-off radius 
2.5a where a is the diameter of the atoms. The param- 
eters e and a are taken to be the same for both phases. 
The parameter c a p — 1 if a = /3; c a p = c otherwise where 
c controls the wetting interaction between the fluid and 
the solid nanoparticle. In this work, we shall consider 
the three cases: c = 0.5, 1, which correspond respectively 
to solvophobic and solvophilic interactions, and c = 2, 
which describes strong bonding between the solid and the 
fluid. The solid nanoparticle is obtained from a spherical 
cut of an equilibrium FCC lattice. In addition to the 
Lennard- Jones interactions, atoms inside the particles 
are connected to their neighbors with FENE springs [?J 
V{r) = -0.5fci? 2 ln(l - (r/R ) 2 ) with k = 30e/cr 2 and 
Rq = 1.5(7. This nearest neighbor bonding allows one 
to heat up the nanoparticle without observing melting 
or fragmentation. This simple modeling is intended to 
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mimic the situation of metallic or oxide nanoparticles, 
with rather high cohesive energies, in volatile organic sol- 
vents with much lower cohesion and a low boiling tem- 
perature. 

Throughout lengths, energies and times are expressed 
in units of a, e and r = ^Jma 2 Je where m denotes the 
common mass of the atoms. For liquid Argon, the cor- 
responding values are: a — 0.3 nm, e = 0.025 eV and 
t — 1 ps. We integrate the equations of motion using a 
velocity Verlet algorithm with a time step dt = 0.005r. 
All the systems considered have been first equilibrated 
at a constant temperature Tq — 0.75 under the constant 
pressure Pq = 0.015 (using a Nose/Hoover temperature 
thermostat and pressure barostat |9j). The temperature 
Tq is below the boiling temperature, that we found to be 
Tfc ~ 0.8, using independent simulations of a liquid/vapor 
interface, under the pressure Pq we are working at. Af- 
ter 100000 time steps of equilibration, the nanoparti- 
cle is heated up at different temperatures T p > Tj, by 
rescaling the velocities of the solid particles at each time 
step, while the whole system is kept at the constant pres- 
sure Pq using a NPH barostat. The fluid beyond a dis- 
tance 10cr from the particle surface is thermostatted at 
T = 0.75, again using velocity rescaling. A global set 
up of the system is depicted in fig. [1] Temperature, den- 
sity and pressure fields have been obtained by averaging 
the corresponding quantities during 10000 time steps in 
nanoparticle centered spherical shells of width ~ 0.15er, 
after a steady state is reached. The Irving-Kirkwood for- 
mula [10l | is used to calculate the normal component of 
the pressure tensor P rr . Finally, we calculate the heat 
flux density flowing through the solid particle, by mea- 
suring the power supply needed to keep the nanoparticle 
at the target temperature T p . 




FIG. 1: Schematic illustration of the simulation setup for the 
nanoparticle and planar cases 

Figure [2] displays steady state temperature profiles 
close to the nanoparticle surface, for different temper- 
atures T p of the nanoparticle. For low T p , the temper- 
ature field in the liquid is practically indistinguishable 
from the form A/r + B, predicted by continuum heat 
transfer equations in homogenous media in spherical ge- 
ometry. Inside the solid, the temperature is not uniform 
but slightly curved downwards, due to the finite conduc- 
tivity of the nanoparticle. Noticeably, the temperature 
is not continuous at the solid interface. The correspond- 
ing temperature jump AT is related to the existence of 
a finite interfacial resistance R = AT/j, j being the flux 
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FIG. 2: Steady state temperature field across the liq- 
uid / nanoparticle interface, obtained with molecular dynam- 
ics. The nanoparticle/fluid interaction is hydrophilic (c = 1). 
The position of the nanoparticle surface is indicated by dashed 
lines, r measuring the distance to the center of the nanopar- 
ticle. Solid curves correspond to fits by the continuum theory 
result T(r) = A/r + B. Note that for T p = 1, the solid curve 
is indistinguishable from the simulation data. 



density flowing through the interface. This effect is well 
known since the pioneering work of Kapitza on Helium, 
and is particularly important when the dimensions of the 
system considered are comparable to the Kapitza length 
Ik = Ai?, A denoting the thermal conductivity of the liq- 
uid 1121] . For usual liquids, Ik is on the order of a few nm 
[ill . Il5l | , thus the effects are particularly important for 
heat transfer around nanoparticles. From fig [21 we mea- 
sured a value R ~ 1.6, which is consistent with previous 
MD simulations @, [H| and consistent with the experi- 
mental determinations [TJ, [l5| in the case of hydrophilic 
interactions, considering that a value of R = 1 in our 
units corresponds to an interfacial resistance on the or- 
der of 0.1 K.m 2 GW _1 . For solvophobic interactions, 
we have measured a larger interfacial resistance (R = 13 
for c=0.5) while for c = 2, no temperature jump is seen 
due to very good thermal contact between the nanopar- 
ticle and the fluid. Upon increasing the temperature of 
the nanoparticle, deviations from the 1/r behaviour are 
clearly seen on fig. [21 when the local temperature exceeds 
T ~ 1, which is about equal to the critical fluid tempera- 
ture. Interestingly, the temperature profile steepens close 
to the nanoparticle surface, corresponding to a decrease 
of the local effective conductivity, and providing "ther- 
mal shield" for the fluid away from the particle surface 
at distances one particle radius. This thermal barrier 
indicates that the particle can be heated to very large 
temperatures, while the liquid at particle diameter away 
from the surface can have close to ambient temperature. 
Also the temperature gradient at the particle surface can 
reach enormous values 0.3 in reduced units (see Fig. 2) 
which translates to 300 K/nm for Argon parameters. To 
relate these features to possible structural changes of the 
fluid, we have considered the density profile around the 
nanoparticle in fig. [3] When heating is moderate, some 
layering is present, reminiscent of the layering observed at 
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FIG. 3: Density profile close to a hot nanoparticle having a 
temperature T v . The inset displays the normal pressure close 
to the hottest nanoparticle (T p = 5). 
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FIG. 4: Flux density as a function of the temperature of 
a nanoparticle (filled symbols) compared to the flux density 
flowing through a flat interface (open symbols). Lines are 
guides to the eye. 



a hydrophilic interface in the absence of heating. When 
the temperature of the particle increases, the height of 
the first peak decreases while layering becomes blurred, 
and a dilute layer appears. Note however that the cor- 
responding density within this layer is still one order of 
magnitude larger than the vapor density at coexistence, 
that we determined using independent simulations of a 
stable liquid/vapor interface. Correspondingly, the inset 
of figure [3] shows a gradual increase of the normal pres- 
sure in the vicinity of the particle, with values of the 
pressure approaching the critical pressure, estimated to 
be P c = 0.1 [HI]. This increase in the local pressure, 
which is similar to a Laplace pressure effect in capillar- 
ity, appears to prevent formation of a vapor. It is also 
worth noting that the decrease of density is accompa- 
nied by the existence of a limiting temperature profile 
far away from the particle, as seen in fig. [21 It is inter- 
esting to contrast this situation with the one obtained in 
a planar geometry. To this end, we performed a set of 
simulations in which a FCC solid slab cut along the 100 
direction is in contact on both sides with the fluid. The 
system is periodic in all directions, and both the solid 
slab and the liquid at a distance 12a from the solid are 
thermostatted at T = 0.75. The system is kept at con- 
stant pressure Pq = 0.015 in the direction normal to the 
solid slab. After equilibration, the temperature of the 
solid slab is raised in order to establish a heat flux be- 
tween the solid and the bulk liquid. The two situations 
(nanoparticle and flat solid) are compared in figure [H in 
terms of the heat flux density as a function of the tem- 
perature of the solid. In each case, the three different 
wetting conditions (c = 0.5, c = 1, c = 2) are consid- 
ered. The difference between the nanoparticle case and 
that of a flat solid is striking. For the nanoparticle, two 
regimes can be distinguished, depending on the heating 
intensity. For low temperatures, the flux increases lin- 
early with T p , up to a kink temperature Tfc. The initial 
slope increases as the wetting becomes better, while the 
kink temperature moves to higher values with increas- 
ing the solvophobic character of the particle. Beyond 



the flux levels off at a plateau value which increases with 
the wetting parameter c. Interestingly, the transition be- 
tween these two behaviors takes place when deviations 
from the 1/r profile become significant, and the liquid 
density of the first peak close to the particle decreases 
significantly. This confirms that the nature of the fluid in 
the interfacial zone is modified significantly for such heat 
fluxes, with a decreased conductivity - or equivalently a 
high value of the effective interfacial resistance. For flat 
interfaces, on the other hand, a quite different scenario 
appears: for low fluxes, the flux increases nearly linearly 
with temperature and the curves are practically indistin- 
guishable from their nanoparticle counterparts. However, 
the heat flux density drops abruptly when the tempera- 
ture is raised above a critical temperature that increases 
with hydrophobicity. The corresponding maximum heat 
flux increases with wetting. This maximum in the heat 
flux can be described as a "critical heat flux" [17| phe- 
nomenon, observed here at the atomic scale. Note also 
that this drop occurs below the kink temperature Tfc, but 
still beyond the boiling temperature X& = 0.8. 

From a microscopic standpoint, the drop in the heat 
flux is associated with the formation of a vapor layer and 
a nonequilibrium drying of the surface. Due to its low 
conductivity, this layer completely blocks heat transfer 
from the solid. The density within the flat vapor layer is 
~ 0.02, one order of magnitude smaller than the dilute 
fluid layer surrounding hot nanoparticles. Furthermore, 
the pressure within this layer is uniform and equal 
to the imposed pressure Pq- This is in contrast with 
the nanoparticle case, where the local curvature of the 
iso-density curves imposes an increase in the pressure, 
and in turn prevents the formation of a well defined 
vapor layer. In spite of this difference, a "critical heat 
flux" can be defined both for the nanoparticle and for 
the planar walls. Any attempt to transfer a flux density 
higher than the limiting values shown in figure [4] would 
result in an unbounded increase of the temperature of 
the nanoparticle, and eventually to its destruction. The 
critical heat flux, however, is increased by a factor of 



4 



almost 4 compared to the planar case. 

Many experiments that probe heat transfer at the 




time (LJ units) 

FIG. 5: Cooling kinetics for a wetting nanoparticle (c = 1), 
after a steady state was reached, where the particle's temper- 
ature was maintained fixed at a temperature T v . 

nanoscale are based on studies of time resolved studies 
of heat transfer in a transient, rather than stationary, 
regime. We therefore consider the influence of the 
critical heat flux on the behaviour observed in such 
experiments. Our simulation of the transient regime 
proceeds by switching off the heating of the nanoparticle 
and monitoring the cooling process into the thermostat- 
ted fluid. Fig. [5] displays the cooling kinetics of particles, 
starting from various initial temperatures T p . Two 
behaviors are to be distinguished. When T p is smaller 
than the kink temperature Tk discussed before (i.e. 
the initial flux is smaller than critical) , the relaxation 
is typically exponential and the interfacial resistance 
extracted from fitting the relaxation curves Q is in 
good agreement with the value found from the steady 
state temperature profiles. On the other hand, when T p 



exceeds Tk, the relaxation is no longer exponential but 
initially linear, which corresponds to a constant heat 
flux flowing through the interface. This early relaxation 
is followed by a second step, where the decay is more ex- 
ponential, with a decay time typically 20% smaller than 
the value reported for lower initial temperature. Hence, 
the change of behavior revealed by our simulations and 
the existence of a critical heat flux around nanoparticles 
could be assessed experimentally by monitoring the 
cooling kinetics in pump-probe experiments. 
In summary, we have investigated heat transfer from 
nanoparticles and planar interfaces using molecular 
simulations, under conditions of high flux and close 
to the boiling transition of the carrier fluid. In both 
cases the existence of a critical heat flux, above which 
heat transfer is impossible, is observed. This heat flux 
is considerably higher in the case of nanoparticles, as 
the formation of a blocking vapor layer is prevented by 
the increased pressure around the nanoparticle. This 
suggests that nanoparticles could be heated to rather 
high temperatures inside the host fluid. The existence of 
the critical heat flux would manifest itself by a change 
of regime in the transient absorption experiments under 
high flux/high temperature conditions Our study 
also suggests that nanoscale features could significantly 
modify heat transfer properties of solid surfaces in this 
regime. 

All the simulations were done us- 
ing the MP code LAMMPS 

( http : / / www. cs . sandria. gov sj plimp/lammps . ht ml ) . 
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